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The incretin effect, which is a unique stimulus of insulin secretion in response to oral 
ingestion of nutrients, is calculated by the difference in insulin secretory responses from an 
oral glucose tolerance test (OGTT) and a corresponding isoglycemic intravenous glucose 
infusion (IIGI) study. The OGTT model of this study, which is individualized by fitting the 
glucose profiles during an OGTT, was developed to predict the glucose profile during an IIGI 
Study in the same subject. Also, the model predicts the insulin and incretin profiles during 
both Studies. The incretin effect, estimated by simulation, was compared with that 
measured by physiologic studies from eight human subjects with normal glucose tolerance, 
and the result exhibited a good correlation (r> 0.8); the incretin effect from the simulation 
was 56.5% ± 10.6% while the one from the measured data was 52.5% ± 19.6%. In 
conclusion, the parameters of the OGTT model have been successfully estimated to predict 
the profiles of both OGTTs and IIGI studies. Therefore, with glucose data from the OGTT 
alone, this model could control and predict the physiologic responses, including insulin 
secretion during OGTTs and IIGI studies, which could eventually eliminate the need for 
complex and cumbersome IIGI studies in incretin research. 

Keywords: Incretins; Oral Glucose Tolerance Test (OGTT); Isoglycemic Intravenous Glucose 
Infusions (IIGI); Individualized Model 

INTRODUCTION 

Incretins are hormones that are secreted from the gastrointestinal tract in response to 
nutrient ingestion; the glucose-dependent insulinotropic polypeptide (GIP) and glu- 
cagon-like peptide- 1 (GLP-1) are the two incretin hormones that have been identified 
(1, 2). The incretin hormones contribute to 50%-70% of insulin secretion in subjects 
with normal glucose tolerance after an oral glucose load (3), while they contribute only 
10%-30% in patients with type 2 diabetes (1, 4). Incretin-based therapy to enhance the 
incretin effect is becoming more and more popular in the treatment of type 2 diabetes 
(5). To calculate the incretin effect in clinical studies, both an oral glucose tolerance 
test (OGTT) and a corresponding isoglycemic intravenous glucose infusion (IIGI) 
study are required. The latter reproduces the blood glucose profiles during an OGTT 
by intravenous glucose infusion (3). However, the IIGI study is quite complex and re- 
quires highly experienced researchers. Therefore, the development of a model that can 
calculate the incretin effect from the OGTT alone will certainly be of significant value. 

Mathematical models of the glucose dynamics have been used to represent the rela- 
tionship mainly between glucose and insulin (6-9). The glucose profiles during OGTTs 
are difficult to predict due to highly variable factors such as the rate of glucose absorp- 
tion from the intestine (10) and the temporal delay of insulin action (11, 12). Further- 
more, the incretin hormones have been recognized as a major insulin secretory stimu- 
lus after an oral glucose load (1,3, 4); and several models that include the incretin con- 
cept have been developed to describe the glucose dynamics during the OGTT in a 
more physiologic manner (10, 13-15). Although those models have some limitations 
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originated mainly from limited clinical data on incretin hor- 
mones, they encompass the relationships among incretin hor- 
mones, insulin, and the rate of the gastric emptying (10, 13, 15). 

In this study, OGTT data from each subject are fitted to an 
individualized OGTT model and the parameters of the individ- 
ualized model are estimated with the initial conditions deter- 
mined by the leave-one-out cross-validation (LOOCV) method. 
The simulation results are validated by the correlation coeffi- 
cients of glucose, insulin and incretin profiles. The ratio of the 
incremental area under the curves (iAUC) are also measured 
between the estimated data from simulation and the measured 
one from clinical studies. From the validation, the proposed mo- 
del in this research is able to represent the physiologic changes 
during OGTTs and IIGI studies, which leads to the successful 
prediction of the incretin effect by simulated OGTTs and IIGI 
studies for the first time. 

MATERIALS AND METHODS 

Study subjects 

Eight healthy subjects with normal glucose tolerance participate 
in the study. None of them are taking medications and none are 
suffering from heart, liver, nor kidney disease or systemic dis- 
eases including diabetes mellitus. The characteristics of each 
subject are shown in Table 1. 

Tests and measurement 

All tests are performed in the morning after an overnight fast 
( > 10 hr). An OGTT is performed by ingestion of 75 g glucose 
solution. With at least 1 week between the procedures, an IIGI 
study is performed in order to duplicate the plasma glucose 
profile during the corresponding OGTT. An intravenous cathe- 
ter is placed in the antecubital vein for blood sampling. For IIGI 
studies, another intravenous catheter is placed in the basilar or 
antecubital vein in the contralateral arm for intravenous glu- 
cose infusion. To measure the plasma glucose concentration, 
approximately 0.5 mL of blood is drawn every 5 min during the 

Table 1 . Characteristics of tlie subjects 



first 2 hr and thereafter every 20 min for the last 1 hr. In addi- 
tion, 9 mL of blood is taken at 0, 15, 30, 60, 90, 120, and 180 min 
for the measurement of plasma insulin and GIF. Plasma glu- 
cose concentrations are determined by the standard glucose 
oxidase method (YSI 2300 STAT plus analyzer, YSI Inc., Yellow 
Springs, OH, USA), and the plasma insulin concentrations are 
measured by a chemiluminescent immunoassay with a Sie- 
mens Centaur XP (Siemens Healthcare Diagnostics Inc., Tarry- 
town, NY, USA). The total GIF concentrations are measured by 
an enzyme-linked immunosorbent assay (Millipore, Billerica, 
MA, USA) based on the manufacturer's instruction. 

Mathematical model and assumptions 

A mathematical model of glucose dynamics encompassing the 
incretin concept has been proposed by Brubaker et al. (13) which 
has been used to simulate the responses to 50 g or 100 g oral 
glucose load. The model consists of three ordinary differential 
equations representing the changes in plasma glucose, insulin, 
and the incretin hormones during the OGTT and other three 
equations representing the absorption of glucose from the gas- 
trointestinal (GI) tract and the glucose balance in the liver. In 
this study, the equations of the GI tract have been modified pro- 
portionally to 75 g glucose load and those are shown in Eq. 1 and 
Eq. 4. 

Duodc (75 g load) = ( 0, t< 5 min 

< 5.1014-0.03071, 5<t< 166.1 min{l) 
i 0, t> 166.1 min 

ll'^'Duodc (75 g load)dt = 414.27 mM = 74.58 g (2) 
RaGutG(75gload)=j 0, t<5min 

1 0.3072a-5/^V^^^^^-^^, t > 5 min 
i^RciGutG (75gload)dt = 413.13 mM = 74.37 g (4) 
The equations describe the delivery of glucose to the duode- 
num (Duodc; mM/min), which in turn determines the rate of 
secretion of incretin hormones and the rate of the absorption of 
glucose into mesenteric circulation (RacutG,* mM/min), which 
determines the changes in blood glucose kinetics. In the IIGI 
study, Eq. 1, 2, 3, and 4 are set to zero. 



Cliaracteristics 



Subjects 



Mean + S.D. 





SI 


S2 


S3 


S4 


S5 


S6 


S7 


S8 




Sex (M/F) 


M 


M 


M 


M 


M 


F 


M 


F 




Age (yr) 


51 


35 


32 


62 


44 


58 


30 


37 


43.6 + 11.4 


HbAlc(%) 


5.5 


6.1 


5.7 


5.7 


5.9 


5.7 


5.3 


5.7 


5.7 + 0.2 


Weight (kg) 


82.2 


74.8 


66.5 


67.2 


49.6 


48.8 


61.5 


55.9 


63.3 + 11.0 


BMI (kg/m^) 


25.7 


24.5 


22 


24.5 


18.8 


20.2 


19.3 


20.2 


21.9 + 2.5 


Incretin effect (%) 


51.2 


49.2 


32.9 


38.5 


86.9 


75.9 


51.4 


33.8 


52.5 + 19.5 


Plasma glucose level of 75 g OGTT (mM/L) 




















Fasting 


4.7 


5.6 


4.9 


5.1 


5.1 


5 


5.1 


5.1 


5.1 + 0.3 


Peak level 


10.2 


8.8 


11 


10.1 


8.4 


9.9 


9.1 


9.7 


9.7 + 0.9 


Peak time (min) 


55 


45 


45 


55 


50 


40 


40 


40 


46.3 + 6.4 


at 120 min 


6.8 


6.4 


7.7 


7.3 


7.3 


7.3 


5.6 


6.7 


6.9 + 0.1 



S.D., standard deviation; M, male; F, female; HbAlc, hemoglobin Ale; BMI, body mass index; OGTT oral glucose tolerance test. 
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The hepatic glucose balance (Hepbalc; mM/min) reflects the 
sum of glucose production in the liver and glucose uptake from 
the mesenteric circulation described in Eq. 5. 

Hepbalc - (HepbalcB, G<Gb 

<HepbalGB+M(GB-G)l G > Gb^I > 0 ^^^^ 

\HepbalGB+M(GB-G)I+a(^- XG> Gb4I < 0 
Gi Gb IB 

HepbalcB (mM/min) is the basal hepatic glucose balance in 
the fasting state (13). M (L^/mU/min) represents the effects of 
counter-regulatory factors on the hepatic glucose balance de- 
termined by parameter estimation in the current study, and a 
(mM^ mU/L^/min) represents the hepatic regulatory term for 
hypoinsulinemia. Gb (mM/L) is the basal glucose level measur- 
ed by the measured fasting plasma glucose level for each sub- 
ject, and Ib (mU/L) is the basal insulin level determined by the 
average fasting plasma insulin level of all subjects. 

The change of plasma glucose concentrations (G; mM/L) is 
described as follows: 

(RacutG Hepbalc i r-13 1 t dl , Giv r-^ m n^r i 

,^ + — ki G-k2 1+y —J- +kio -r— , G< 10 mML 

dG_J V V at V 

1^ G"-fe/+y ^- ^^ + fc»-^, G > lOmML' 

^ V V at V V 

Where V (L) is twenty percent of the weight of each subject, 
ki (L°^/mM° ^/min) and k2 (mM/min/mU) represent non-insu- 
lin-mediated glucose uptake and insulin-mediated glucose up- 
take, respectively. Parameter y is a shaping factor for the deriva- 
tive control of insulin on glucose, ks (L/min) and Iq (mM/min) 
are the slope and intercept of renal glucose clearance, respec- 
tively, and Giv (mL/min) represents the intravenous glucose in- 
fusion rate of the IIGI studies. The term Giv has been added to 
represent the intravenous infused glucose, a term that is used 
in the simulation for the IIGI study, but not for the OGTT. Pa- 
rameter kio (mM/mL) is the correction factor because the unit 
of the intravenous glucose infusion rate are mL/min (glucose: 
180 g/M, 1.54 g/cm^), and the initial value is set to 10. 

The incretin concentration (Inc; ng/L) is described in Eq. 7 on 
the basis of the GIP concentrations during the OGTT. 

dMc ^ Ra^ +k5DuodG-keInc (7) 
dt V 

Rainc = keVIncB (8) 

Rainc (ng/min) is the basal rate of appearance of the incretin, 
IncB (ng/L) is the basal incretin concentration determined by 
the average fasting total GIP of all study subjects and parameter 
ke (min^) is the degradation of Inc. Parameter ks (ng/L/mM) is 
the rate of appearance of incretin due to DuodG, which is de- 
termined by parameter estimation in the current study. In the 
IIGI studies, the rate of appearance of Inc is nearly unchanged 
because the rate is most affected by Duodc. 

The insulin concentration (I; mU/L) is described in Eq. 9 to 
reflect the effect of both glucose and incretin on the insulin se- 
cretion. 



^kyG^ '+ksInc-ksI+fi (9) 

Parameters k? (mU/min/mM^^/L^^) and ks (mU/min/ng) are 
the rate of appearance due to glucose and incretin, respectively 
and those are determined by parameter estimation in this study. 
Parameter kg (min^) is the clearance rate of insulin, and (3 (mU/ 
L/min) is the effect of additional regulators. In the IIGI study, 
the insulin secretion by incretin is simulated to be much lower 
than that of simulated OGTT. 

The model is implemented by Matlab R2012a (Mathworks 
Inc., Natick, MA, USA) along with parameter estimation. The 
solutions to the differential equations in the model are obtain- 
ed by the 'ode45' function, which is based on the 4-5th order 
Runge-Kutta method. The simulation for OGTT shows the re- 
sponses depending on the oral route of glucose administration 
by Eq. 1 and Eq. 3. In contrast, the simulation for the IIGI stud- 
ies shows the responses following the intravenous route of glu- 
cose administration by the term related to GIV in Eq. 6. The in- 
travenous glucose infusion rate in the simulation is determined 
using the same interval as the IIGI study without any delay. The 
initial conditions for glucose (G), insulin (I), and incretin (Inc) 
for both simulations are obtained from the OGTT data. 

Parameter estimation 

The objective of parameter estimation is to fit the profiles of mea- 
sured data. Following the process described in Fig. 1, the values 
of the parameters of the OGTT model are estimated, and the 
probability distribution of each parameter based on the mea- 
sured data is determined. To obtain the appropriate range of 
the parameters, an LOOCV method is utilized prior to estimat- 
ing the parameters of the individualized model. The LOOCV 
method is suitable for the analysis of a small data set and it has 
been shown that LOOCV provides a fairly unbiased estimation 
of the true generalization ability of a classifier (16). In this study, 
the measured data from OGTT are divided into single versus 
remaining classes. A single observation means the OGTT data 
of one subject, which is used to estimate the values of the pa- 
rameters of an individualized OGTT model to fit the actual OG- 
TT data of a subject. The remaining observations mean the OG- 
TT data of all the other subjects except one, which is used to set 
a probability distribution and the initial conditions of the pa- 
rameter estimation for the individualized OGTT model. This 
process is repeated for each study subject. 

The values of the parameters of the OGTT model are estimat- 
ed to minimize the weighted sum of the squares of the residuals 
between the measured and simulated results of the OGTT (17). 
Depending on the objective, different types of residuals are se- 
lected. The plasma glucose (G), insulin (I), and GIP (Inc) data 
are used to set the probability distribution using the remaining 
observations, whereas only the plasma glucose (G) level is used 
to estimate the values of the parameter of the OGTT model. 
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The four parameters ks, k?, ks, and M among the fifteen pa- 
rameters in the mathematical model are estimated with directly 
measured data. The estimation has a limit within a reasonable 
range with the nonzero glucose and the insulin levels and cov- 
ers the parameter distribution. The initial conditions for glucose 
(G), insulin (I), and incretin (Inc) for parameter estimation are 
obtained from the OGTT data, and the initial values for the four 
parameters are obtained from the result of the LOOCV method. 

At first, parameter kio is set to a fixed value, but after estimat- 
ing other parameters, the value of kio is estimated to minimize 
the sum of squares of the differences between the measured 
and simulated results of glucose response during the IIGI study. 
The range of the value is between 6 mM/mL and 14 mM/mL. 

The results of parameter estimation are obtained by the 'pat- 
ternsearch' function (Matlab built-in function) which is based 
on the pattern search method for optimization techniques. 

Validation 

In this research the predicted result with simulation is validated 



All subjects 
(n = 8) 



75 g OGTT 



Repeat n times 



IIGI study 



Measured 
data 



Leave one 
model 



out for 



n-1 

measured data 



1 

measured data 



Parameter estimation 
with G, I, and GIP data 



Initial condition, 
range of parameters 



V 



Parameter estimation 
with G data 



OGTT model 



Intravenous glucose 
infusion rates for subject 



Simulation of 




OGTT & IIGI 


Simulation 



results 



Validation 



1) Correlation coefficient 

2) Ratio of lAUC 

3) Incretin effect 

4) Root mean square error 

Fig. 1. A schematic overview of the process of parameter estimation for the model 
and validation between measured data and simulation results. OGTT, oral glucose 
tolerance test; IIGI, isoglycemic intravenous glucose infusions; G, glucose; I, insulin; 
GIP, glucose-dependent insulinotropic polypeptide; lAUC, incremental area under the 
curves. 



in several ways. First, Pearson's correlation coefficient between 
measured and simulated data is used as a measure for the good- 
ness of fit and the accuracy is evaluated with the ratio of the iAUC 
of each profile between simulated and measured results (ratio = 

iAUC of the profile of simulation 
iAUC of the profile of measured data 

Second, the incretin effect is calculated from insulin secre- 
tion in the simulation result and then it is compared with the 
value from the measured data (1, 2). It is calculated by the dif- 
ference of the iAUC between the insulin profile of an OGTT and 
its corresponding IIGI study by the trapezoidal rule - 

iAUCofinsulinprofflesofOGTr-iAUCofinsulinprofflesofllGIstudy ,^ 

(Incretin effect (%) = 100 x J. 

iAUC of insulin profiles of OGTT 

Lastly, the intravenous infusion rates used during the IIGI 
study for each subject are applied to the corresponding OGTT 
model to validate the results of parameter estimation. A high 
correlation coefficient and a low root mean square error (RMSE) 
represent the most suitable results. For example, if the intrave- 
nous infusion rates used for subject #1 are applied to the eight 
individualized models, the IIGI simulation with parameters 
from the OGTT model for subject #1 should have the highest 
correlation coefficient and the lowest RMSE compared with the 
actual IIGI study data of subject #1. 

Ethics statement 

The research protocol has been approved by the institutional 
review board of Seoul National University Hospital (IRB No. 1011- 
086-341). All participants were informed of the experimental 
procedures and gave their informed consent before any study- 
related procedures. 

RESULTS 

The initial conditions and the range of estimated parameters 
were determined to reflect the characteristics of measured data. 
The initial conditions and probability distribution determined 
by the remaining measured data described in Fig. 1 had an im- 
pact on the results of the parameter estimation. The most no- 
ticeable discrepancy in the measured data between the data in 
this study and the data reported by Brubaker et al. (13) was an 
attenuated secretion of insulin and GIP. This characteristic was 
reflected in the results of estimation as ks and k? in Table 2. Pa- 



Table 2. Parameter estimation results of tlie OGTT model 



Parameters 


Description 


Source 


Estimated value* 


Reported value (1 3) 


Units 


ks 


Rate of appearance of incretin due to Duodc 


Derived from steady state 


6.91 + 0.0029 


27.64 


ng/L/mM 


kz 


Rate of appearance of insulin due to G 


Adjusted parameter 


0.095 + 0.024 


0.125 


mU/min/mM^7L°-2 


ks 


Rate of appearance of insulin due to Inc 


Adjusted parameter 


0.0042 + 0.0016 


0.0050 


mU/min/ng 


M 


Effects of counter-regulatory factors on the liver 


Adjusted parameter 


0.030 + 0.009 


0.020 


L7mU/min 



^Estimated values are presented as mean : 
tration; Inc, plasma incretin concentration. 



standard deviation. OGTT, oral glucose tolerance test; Duode, the delivery of glucose to the duodenum; G, plasma glucose concen- 
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Table 3. The validation results between simulation and measured data 



Variables 




OGTT 






IIGI study 




Glucose 


Insulin 


Incretin 


Glucose 


Insulin 


Incretin 


Correlation coefficient (mean + S.D.) 


0.81 + 0.17 


0.85 + 0.18 


0.59 + 0.42 


0.85 + 0.06 


0.85 + 0.11 


0.18 + 0.56 


Correlation coefficient for #7 


0.91 


0.87 


0.13 


0.96 


0.93 


0.02 


Ratio of lAUC* (mean + S.D.) 


1.18 + 0.39 


0.79 + 0.26 


0.86 + 0.41 


1.14 + 0.35 


0.86 + 0.52 


0.67 + 0.86 


Ratio of lAUC for #7 


1.25 


0.93 


1.21 


1.06 


0.93 


0.24 



*Ratio of lAUC, the ratio of the incremental area under the curve for profiles between the simulation results and the measured data. OGTT oral glucose tolerance test; 
glycemic intravenous glucose infusions. 



:il, iso- 



15 



1,000 




0 30 60 90 120 150 180 



Time (min) 

Fig. 2. Plasma glucose concentration profiles in subject #7. The solid lines are from 
the simulation (blue, 75 g OGTT red, IIGI study; green, 4 glucose infusion rates). The 
dotted lines are from the measured data (sky blue [•], OGTT purple [o], IIGI study). 



60 r 



50 - 



^ 40 - 




Time (min) 

Fig. 3. Plasma insulin concentration profiles in subject #7. The solid lines are from 
the simulation (blue, 75 g OGTT red, IIGI study). The dotted lines are from the mea- 
sured data (sky blue [•], OGTT purple [o], IIGI study). 

rameter ks, describing the rate of appearance of incretin hor- 
mones due to Duodc, and parameter k?, describing the rate of 
appearance of insulin due to G, were lower by 75% and 25%, re- 
spectively, compared with the value (13). Parameter M, which 
was affected by the hepatic glucose balance, and parameter ks, 
which described the rate of appearance of insulin due to incre- 



800 - 



f 600 - 
o 




0 30 60 90 120 150 180 



Time (min) 

Fig. 4. Plasma incretin (GIP) concentration profiles in subject #7. The solid lines are 
from the simulation (blue, 75 g OGTT; red, IIGI study). The dotted lines are from the 
measured data (sky blue [•], OGTT; purple [o], IIGI study). 

tin, were 30% and 16% higher compared with the value (13), re- 
spectively (Table 2). 

The correlation coefficient and the ratio of iAUC between the 
simulation result and the measured data of insulin were over 
0.80 (Table 3). The correlation coefficient and the ratio of iAUC 
are also shown in Table 3. The mean value of correlation coeffi- 
cient of incretin was smaller than those of insulin and glucose. 
The mean ratio of iAUC of glucose was greater than 1, while those 
of insulin and incretin results were less than 1. For a representa- 
tive example, the profiles for the simulation results of subject #7 
from the OGTT and IIGI study and the directly measured data 
(with symbols) are depicted in Fig. 2, 3, 4 and the correlation 
coefficient and the ratio of iAUC for subject #7 are shown in Ta- 
ble 3. 

For the first simulation of the IIGI study, parameter kio was 
fixed at a value of 10 for all subjects. However, the glucose ex- 
cursion of each subject was varied during the IIGI study. Among 
the various conditions that influence the glucose homeostasis, 
delayed insulin response to the glucose administration was fo- 
cused. The response of the intravenously infused glucose was 
increased by 15.44 ± 9.64 mg/dL per 1 mL/min in average at 
the initial time of the IIGI study because of the time lag of insu- 
lin secretion. Despite these variations, the response to sequen- 
tial intravenous glucose infusion could be predicted by using 
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the mathematical model. For the second simulation, kio was es- 
timated to minimize the difference between the measured data 
and the simulation, reflecting the different responses to intra- 
venous glucose administration for each subject. The estimated 
value of kio (mean ± S.D.) was 9.25 mM/mL ± 3.06 mM/mL. The 
range of the value was from 6 mM/mL to 13.3 mM/mL. The cor- 
relation coefficient of the glucose profiles using the estimated 
kio value for the IIGI studies was higher at 0.87 ± 0.08 than that 
using the fixed value. The RMSE of the glucose profiles using the 
estimated fco value was reduced by 50% to 360.6 mM/L ± 440.3 
mM/L compared with a value of 721.9 mM/L ± 815.0 mM/L with 
the fixed value. 

DISCUSSION 

In this study, measured data of the subjects have a difference 
with other reported data used OGTT model. The difference was 
the low secretions of insulin and GIP which might be caused by 
different experiment and different characteristics of study sub- 
jects (2, 4, 6, 12, 18-23). Due to the heterogeneity among human 
studies, the range of parameters can hardly be determined on 
the basis of earlier studies. So, LOOCV method described in me- 
thod has been applied to estimate initial conditions of parame- 
ters in order to reflect the difference in the measured data be- 
tween the current study and a report by Brubaker et al. (13). The 
estimated initial conditions and the range of parameters helped 
the OGTT model to closely fit the individual profiles. 

The lower secreted insulin does not affect the 'incretin effect' 
in the subject. Despite the apparently lower insulin responses 
during OGTTs, the incretin effect calculated by both measured 
and simulated data was within the previously reported range, 
i.e., 50%-70% (1, 4, 24): the incretin effect from the simulated 
data was 56.5% ± 10.6% and that from the actual experiment 
was 52.5% ± 19.6%. 

The glucose profiles generated during IIGI simulation based 
on the OGTT model should be examined by another validation 
for suitability. These validation results were shown in Fig. 5. The 
correlation coefficient between glucose profiles during actual 
IIGI and simulated IIGI for subject #7 had a highest value of 
0.96 when the simulated IIGI was adapted from its own OGTT 
profiles. The RMSE for subject #7 was 5.13 mM/L and the val- 
ues for the other subjects ranged from 4.61 mM/L to 17.78 mM/ 
L. The value from subject #7, 5.13 mM/L, ranked the second 
lowest, but this was comparable to the lowest value from sub- 
ject #8, 4.61 mM/L. The reason for this similarity could be ex- 
plained by the fact that these two subjects had similar glucose 
profiles during OGTTs, fasting glucose, time of peak glucose 
level, and glucose response at 180 min as shown in Table 1. The 
other cases also had similar tendencies. Thus, the IIGI simula- 
tion based on the OGTT model could accurately predict the re- 
sponse to intravenous glucose administration. 
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Fig. 5. The IIGI simulation for all OGTT models with the same intravenous glucose in- 
fusion rates for subject #7. The solid lines with different colors show the OGTT model of 
eight subjects (si -8). The black dashed line shows the measured data for subject #7. 

The proposed model have some limitations. Firstly, the ordi- 
nary differential equations of the OGTT model represented only 
general single peak of glucose level during the OGTT study. In 
case of subject #5, unusual double glucose peak levels were re- 
presented during the OGTT, which was not reproducible with 
our model so that the correlation coefficient for glucose profiles 
from subject #5 was very low (0.40). Nevertheless, the IIGI re- 
sponse was well predicted with a correlation coefficient of 0.90 
because the simulation reflected the glucose levels from the in- 
travenous glucose infusion. Secondly, since a limited number 
of subjects were used in this study, the proposed model needs 
to be further validated with additional normal subjects which is 
being planned as a future study. Thirdly, since all the subjects 
in this research have normal glucose tolerance, their responses 
are in normal range. It is necessary to develop advanced model 
based on this study that may predict the incretin effect and in- 
cretin response in the different populations such as subjects with 
type 2 diabetes or abnormal glucose homeostasis. Fourthly, only 
GIP results in simulation were compared to those in clinical 
study even though both GIP and GLP-1 are incretin hormones 
which contribute to the incretin effect (25). However, it was as- 
sumed that the changes in GLP-1 are too small to influence the 
model because plasma levels of GLP-1 are much lower than 
those of GIP in humans (13). Further studies to improve the 
model considering both incretin hormones (i.e., GIP and GLP- 
1) would be necessary. Lastly, fasting glucose levels at 0 min 
and decreased glucose levels near 180 min were unpredictably 
lower or higher than the expected range of parameters. To over- 
come these discrepancies, most of the glucose levels in the OG- 
TT model was adjusted. Moreover, the initial conditions in this 
study were obtained by either OGTT or LOOCV method and 
only the glucose profiles during OGTTs for parameter estima- 
tion was used, although the data for insulin and GIP were also 
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obtained. In this research, the objective was not fitting the mod- 
el to the measured data, but predicting the other profiles using 
the OGTT model. The predicted insulin and incretin profiles 
were actually similar to the measured data, though the initial 
conditions were unavailable prior to the modeling. 

In summary, for the calculation of incretin effect from the dif- 
ference in insulin secretion, the proposed model in this research 
has attempted to represent the different routes of glucose ad- 
ministration; oral vs intravenous. The parameter estimation by 
applying the LOOCV method was able to reflect characteristics 
of the measured data without complicated comparisons. More- 
over, the insulin and incretin profiles were effectively predicted 
to be based on the glucose response depending on the route of 
administration of glucose. The outcomes of this research will be 
useful to control and predict the responses of OGTTs and IIGI 
studies. This is the first proposed model to predict the hormone 
response of IIGI study fi:om the data of OGTT study and the mo- 
del would be the first step to understanding the incretin physi- 
ology with simulation study. Furthermore, by taking advantage 
of the prediction of the proposed model, the accuracy of the IIGI 
study could be improved by predicting the intravenous glucose 
infusion rates to trace the glucose profile during the OGTT. As a 
result, in clinical practice, the proposed method could signifi- 
cantly reduce the burden on subjects and researchers to predict 
the incretin effect. As a conclusion, the current model could 
play a vital role in developing a better understanding of the in- 
cretin effect and the proposed method could be used in addi- 
tional clinical studies for the investigation of glucose dynamics. 
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